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f ^ ' We present a detailed analysis of high resolution observations of the three lowest CO transitions in two nearby active galaxies, 

NGC 4569 and NGC 4826. The CO(1-0) and (2-1) lines were observed with the Plateau de Bure Interferometer and the CO(3-2) line 
with the Submillimeter Array. Combining these data allows us to compare the emission in the three lines and to map the line ratios, 
|j R21 = /co(2-i)/^co(i-0) and = /co(3-2)/^co(i-0) at a resolution of ~2", i.e. a linear resolution of 160 pc for NGC 4569 and 40 pc for 

NGC 4826. In both galaxies the emission in the three lines is similarly distributed spatially and in velocity, and CO is less excited 
(i?32 < 0.6) than in the Galactic Center or the centers of other active galaxies studied so far. According to a pseudo-LTE model the 
molecular gas in NGC 4569 is cold and mainly optically thick in the CO(1-0) and (2-1) lines; less than 50% of the gas is optically 
thin in the CO(3-2) line. LVG modeling suggests the presence of an elongated ring of cold and dense gas coinciding with the inner 
Lindblad resonance (ILR) of the stellar bar in agreement with a previous analysis of the kinematics. More excited gas is resolved in 
the circumnuclear disk of NGC 4826. According to our pseudo-LTE model this corresponds to warmer gas with a ~50% of the CO(3- 
2) emission being optically thin. LVG modeling indicates the presence of a semicircular arc of dense and cold gas centered on the 
dynamical center and ~70pc in radius. The gas temperature increases and its density decreases toward the center. A near side/far side 
asymmetry noticeable in the CO, R32 and Paa maps suggests that opacity effects play a role. Examining published CO maps of nearby 
I active galaxies we find similar asymmetries suggesting that this could be a common phenomenon in active galaxies. These mainly 

, qualitative results open new perspectives for the study of active galaxies with the future Atacama Large Millimeter/submillimeter 

C\) ' ^ray- 
Key words, galaxies: individual: NGC 4569, NGC 4826 - galaxies: active 

T"! | Introduction model for active galactic nuclei (AGN; lAntonuccilll993l) . As 

>- , . . . , . , . , the most abundant molecule after H2 and with its four lowest ro- 

' Multltransltlon CO observations are essential for studying the transitions coinciding with the atmospheric windows in 

^ . processes at work m the centers of active galaxies. The molecu- ^ minimeter and submi i Umeter range , CO has long been rec- 
H ; lar gas fuels the nuclear activity powered by accretion onto a su- ized and used ag ^ begt molecular fracer observing 

permassive black hole and/or intense star formation. Conversely, a ^ transition allows us to the molecular gas col . 

the molecular gas properties are affected by the influence of the umn densi and kinematics ancL thus ^ to identif dynamical in _ 

active nucleus through jets winds, and radiation. The molecular stabiM and measure gravitationa i torques and possible gas 

gas must also constitute the bulk of the nuclear torus responsi- inflows Such deta0ed smdies ^ for e k conducted for 
ble for the highly anisotro pic extinction invoked by the unified by in {h context of ^ NTJdei of GAlaxies survey 

Send offprint requests to: F. Boone (|Garcia-Burillo^taL| |200j hereafter NUGA-I) based on high 

* Based on observations carried out with the Submillimeter Array resolution observations of the CO(l-O) and (2-1) lines (see also 

(SMA), a joint project between the Smithsonian Astrophysical |Helfer et al.||2003t |Sakamoto et al.||1999t for other high resolu- 

Observatory and the Academia Sinica Institute of Astronomy and tion CO surveys). It is however necessary to observe more tran- 

Astrophysics that is funded by the Smithsonian Institution and the sitions to characterize the physical conditions of the molecular 

Academia Sinica. gas and identify interstellar medium (ISM) structures involved 

** Based on observations carried out with the IRAM Plateau de Bure in star formation, shocks (e.g., due to the disk kinematics, AGN 

Interferometer. IRAM is supported by INSU/CNRS (France), MPG feedback, stellar winds, or supernova remnants), or quiescent gas 
(Germany) and IGN (Spain). 
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reservoirs ready to ignite star formation or feed the central black 
hole. Other molecular or atomic lines (including those of CO 
isotopologues) can be observed to characterize the physical con- 
ditions, but this implies adding new unknowns, namely the rela- 
tive abundances of the observed molecules or atoms (which are 
otherwise interesting to measure to study the ISM chemistry). 

Models of the ISM in galactic centers suggest that high res- 
olution multitransition CO observations should allow for the 
identification of the toru s around the supermassive black hole 
dWada & Tomisakal2005l) as well as the influence of the AGN on 
the ISM chemistry (in X-ray dom inated regions, Mei ierink et al.l 
120071 iMeiierink & Spaansll2005l) . 

While the first two CO transitions at 1 15 and 230 GHz have 
been regularly observed for some time at high spatial resolution 
with millimeter interferometers, the third transition at 345 GHz 
could only be observed with single dish telescopes until re- 
cently. CO(3-2) single dish observations complemented lower 
transitions' observations and yielded a wealth of information 
on th e excitation of the molecular ga s in galaxy centers (see 
e.g. Parjadorjoulos & Seaquist 1999; Paoadorjoulos & Allen 
2000; bumke et al.ll200U ICurran et all 120011: lHafok & Stutzkil 
20031: lYao et all 120031: fWeiB et al] 120051; iBavet et al.l 12 006; 
Muraoka et al.l 120071: iKohno et al.l 120081: Papadonoulos et afl 
2008: 1 Wilson et al Il2009t Ilsraelll2009h including that of our own 
Galaxy (lOka et al.ll2007l) . These observations showed that the 
CO(3-2)/(l-0) line ratio generally increases toward the cen- 
ter, where it can reach values greater than unity (in temperature 
units) in starburst and AGN hosts as well as in the Milky Way 
(vs. a typical value of 0.4 in galaxy disk s). It is only sin ce the 
Submillimeter Array became operational dHo et alj|20 04) that it 
has been possible to map the CO (3-2) (and higher /-transitions 
see, e.g., Matsushit a et al . 2009) in the central regions of galax- 
ies with the same spatial resolution as the CO(1-0) and CO(2-l) 
lines observed with millimeter interferometers. Among other re- 
sults, these observatio ns have revealed gas int eracting with a jet 
in the center of M51 dMatsushita et al.l [20041) and gas interact- 
ing with a starb urst superwind or AGN outflow in NGC6240 
(Iono etal . 20071 

In this paper, we present a detailed analysis of high reso- 
lution observations of the three lowest CO transitions in two 
nearby active galaxies. Long-term drivers for this work are 
(1) to estimate the feasibility of mapping the physical con- 
ditions in the centers of galaxies based on CO line observa- 
tions, and (2) to discuss the implications for the study of ac- 
tive galaxies with current and future instruments. We selected 
the two galaxies in the NUGA sample with the strongest low-/ 
CO lines, NGC4569 and NGC4826, and observed the CO(3- 
2) fine with the Submillimeter Array (SMA). The CO(1-0) 
and CO(2-l) observations obtained with the IRAM Plateau 
de Bure I nterferometer (PdBI) were presented and analyzed in 
NUGA-HlGarcfa-Burillo etaTI d2005l hereafter NUGA-IV) and 
Boone et al.l d2007l hereafterNUGA-VID. 

NGC 4569 is a bright early type spiral galaxy in the Virgo 
cluster at a distance of 17 Mpc (1" i s ~80pc). It harbours a 
transition-type nucleus dHo et al . 1997) which exhibits the most 
pronounced nuclear starburst activity among the low-ionization 
nuclear emission-li ne region (LINER) and transition nucle i with 
available UV data dMaoz et al.ll 19981) . iHensler et alJ (120031) re- 
port a bipolar outflow in X-ray an d a giant Ha outflo w with a 
scale height of 10 kpc to the west. IChvzv et al.l d2006l) observe 
symmetric radio lobes extending up to 24 kpc from the galactic 
disk. High resolution PdBI ob servations all owed us to model the 
kinematics of the central kpc dNUGA-VIJ) and showed that the 
molecular gas concentration at a radius of 600 pc corresponds to 
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Fig. 3. Composite color image of NGC 4569 made from the in- 
tensity maps of the CO(1-0), CO(2-l) and CO(3-2) line emis- 
sion. The line intensities are represented by the red, green and 
blue color intensities respectively. The intensity maps were ob- 
tained by summing all channels where the line has S/N>3, and 
are corrected for primary beam attenuation. 



the inner Lindblad resonance (ILR ) of the large scale b ar, which 
extends to a radius of 5 kpc (lLauri kainen & Salo 2002). The in- 
clination of th e galaxy is 70 deg dTullvl Il988h and its position 
angle 30deg (Guhathakur ta et al.l 1988). 

NGC 4826, also known as the 'Evil Eye', is the closest tar- 
get (4.1 Mpc, i.e. 1" is -20 pc) of the NUGA core sample for 
which we have acquired ~0.5" resolution PdBI CO observations 
(iNUGA-IVl) . These observations show a large concentration of 
molecular gas within a radius of 80 pc forming a circumnuclear 
molecular disk (CND). A detailed analysis of the kinematics, 
howe ver, does n ot reveal any evidence of fuelling of the nu- 
cleus dNUGA-j) . The nucleus of NGC 4826 is also classified as 
a transi tion type (|Ho et al.l 19971) . The inclination of th e galaxy is 
60 deg dRubinll 19941) . and its position angle 1 12 deg dNUGA-J) . 

The paper is organized as follows. The SMA CO(3-2) ob- 
servations are presented in Section|2] The distributions of the 
emission in the three lines and of the two line ratios (/?2i = 
/co(2-i)//co(i-o) and R 32 = fco(3-2)/-fco(i-0), where / co is a 
line intensity in temperature units) are compared in Section|3] 
and Section|U respectively. Mapping of the physical condi- 
tions based on radiative transfer modeling of the line ratios is 
performed in Section[5] The results obtained are discussed in 
Section|6]and conclusions are drawn in Section|7] 

2. Observations and data processing 

Each galaxy was observed with the SMA in an ~8 hour track 
at night in March 2005. The observing frequency was set to 
the CO(3-2) line frequency (345.796 GHz) redshifted according 
to the systemic velocity of the galaxy. At these frequencies the 
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Fig. 1. CO(3-2) channel maps of the galaxy NGC4569. The velocity relative to the systemic velocity of the galaxy (vh e i=- 
235kms~') is given at the top left of each map. The phase center indicated by the large cross is at the AGN position 
a/2000=12 /, 36'"49.8 s , (5/2000=13°09'46.3". The beam (represented at the bottom left of each map) is 2.6"x2.1" and the rms is 
115mJybeam _1 . The only contour line overlaid corresponds to the Act level. The maps are not corrected for the primary beam 
attenuation. 



half power beam width (HPBW) of the primary beam is ~36". 
The SMA was in its compact configuration with seven working 
antennas. Saturn was observed for bandpass calibration, bright 
nearby quasars (3C273 and 3C279 for NGC4569, 3C279 and 
1 159+292 for NGC 4826) were observed every ~15 min for gain 
calibration, and Titan and Callisto were used for flux cali bration . 
The data were reduced using the SMA MIR package (Oi 2005). 



For i mage reconstruction the MIRIAD softwar^U (Saul t et al.l 
1 19951) was used with natural weighting applied to the uv data. 

We estimated the flux missed by the SMA due to the short 
spacings inherently missing in interferometric ob servations by 
comp aring to the single dish data publish ed by IWilson et al.l 
(2009) for NGC 4569 and bv llsraell (l2009h for NGC 4826. We 



http://www.cfa.harvard.edu/sma/miriad/manuals/manuals.html 
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Fig. 2. SMA CO(3-2) integrated map of NGC4569 (grey scale given in the color bar in Jy kms -1 beam" 1 ) overlaid with the PdBI 
CO(l-O) integrated map (left) with contours at 2, 6, 12, 20, and 30 Jy kms -1 beam -1 and the PdBI CO(2-l) integrated map (right) 
with contours at at 7, 21, 42, 70 and 87 Jy kms -1 beam -1 . No primary beam correction has been applied. All maps have the same 
resolution (PdBI data have been tapered) and the beam (2.3"x2.1") is shown at the bottom left. The black cross labelled E shows 
the position of the AGN (aj2ooo = 12 h 36 m 49*80, tfnooo = 13°09'46'.'3). It is also the dynamical center and the phase center for the 
observations. The other crosses indicate the positions of the spectra plotted in Fig.|4] 



compared the total flux within an aperture of ~36" in diameter 
after correcting the SMA data for the primary beam attenuation 
and degrading to the resolution of the single dish data (14.5")- 
We found that 21% and 38% of the flux is missed by the SMA 
in NGC 4569 and NGC 4826, respectively. 

We then used the single dish data to compute and add the 
missing short spacings to the SMA data. The corresponding 
large scale distribution was included in the image plane using 
the IMMERGE task of the MIRIAD software. One of the input 
parameters to this task is a factor to scale the flux of the single 
dish data. We controlled the value of this parameter by compar- 
ing the flux of the merged data with that of the single dish data 
following the same procedure as the one described above to es- 
timate the missing flux. The factor was fixed such that the fluxes 
of the merged data and single dish data agreed to within 10%. 

For NGC 4569 the synthesized beam HPBW is 2.6"x2.1" 
with a position angle of -36 deg and the RMS is 1 15 mJy beam -1 
in lOkms -1 channels. For NGC 4826 the HPBW is 2.6"xl.9" 
with a position angle of -37 deg and the RMS is 75 mJy beam -1 
in lOkms -1 channels. 

To make the comparison between the three CO lines easier 
the PdBI observations (which also include the short spacings) 
have been degraded to the resolution of the SMA observations. 
This was done by apodizing the visibilities, cleaning the dat- 
acubes and restoring them with the same clean beam as the one 
used for the corresponding SMA map. 

3. Comparison of the three CO line intensity 
distributions 

3.1. NGC 4569 

The channel maps of the CO(3-2) line emitted by NGC 4569 
and observed with the SMA are presented in Fig.Q] The elon- 
gated distribution identified in most of the individual CO(1-0) 
and CO(2-l) channel maps ( see INUGA- VIII) can also be recog- 
nized in the CO(3-2) channel maps. The CO(3-2) emission ap- 
pears more clumpy and irregular, however, and while emission 



was detected up to +200 km s relative to systemic in CO(l- 
0) and CO(2-l), there is no CO(3-2) emission at a level higher 
than 3 <x for velocities greater than 130 km s - 1 in absolute value. 
These differences between the PdBI and SMA data may be due 
to the lower S/N ratio of the CO(3-2) SMA data. Indeed, the in- 
tegrated map (see RGB image in Fig. [3] and detailed grey scale 
maps in Fig.[2j shows that the CO(3-2) is well detected at the 
peaks of CO(1-0) and CO(2-l) line emission, i.e. the northern 
(B) and the southern (G) peaks that were interpreted as emission 
from an elo ngated ring at the Inner Lindblad Resonance (see 
iNUGA-VIll) . and at the central ridge (D). As in the two other 
line maps, there is a depression in the CO(3-2) map at the loca- 
tion of the nucleus (E) which is also the dynamical center, and 
the southern peak is more prominent than the northern one. 

There are, however, differences in the integrated maps as 
well. The depression at the center and the difference between 
the northern and southern peaks seem to be more pronounced in 
the CO(3-2) map. The CO(3-2) map also shows more emission 
west of the center (F) than the other line maps. 

The spectra at the 8 positions marked by crosses in the inte- 
grated maps (Pig.[2j and labelled from A to H are presented in 
Fig-E] The SMA spectra show negative values outside the CO(3- 
2) line. These are due to deconvolution errors (a single SMA 
track was obtained, whereas the PdBI observations include four 
tracks in four configurations, leading to better uv-coverage and 
therefore lower sidelobes). These effects are difficult to correct 
for without additional observations, but it can be seen that the 
negative offsets do not exceed 100 mJy, i.e., the deconvolution 
errors do not exceed the RMS. 

Although the CO(3-2) line has a lower S/N, the three CO 
lines look fairly similar in their profiles at a given position as 
well as in their spatial distribution, i.e., their intensity ratios are 
constant as a function of velocity to within -30% over the pro- 
files. This is particularly true at the three peaks (positions B, D 
and G). This implies that overall, the kinematics are the same for 
the three emission lines, i.e. even if different clumps or gas com- 
ponents are probed by the different lines they all seem to follow 
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Fig. 4. Individual spectra of CO(l-O) (black), (2-1) (blue) and (3-2) (red) lines at six different positions (see Fig. [2]) in the cen- 
tral region of NGC4569 extracted at the same spatial resolution of 2.6"x2.1". The fluxes have been corrected for primary beam 
attenuation and the CO(l-O) flux is scaled up by a factor 3.5 for better readability. 



the same kinematics. The CO(l-O) emission seems to be rela- 3.2. NGC4826 
tively stronger away from the peaks (positions A, C and H); this 

is particularly striking to the north (at position A), but a similar The channel maps of the CO(3-2) emission from NGC 4826 as 
trend is seen toward the south in position H. observed with the SMA are presented in Fig.|5] The phase cen- 

ter is offset by ~10" to the south of the dynamical center, which 
is marked by a small cross on the channel maps. The typical 
disk rotation pattern (a.k.a. spider diagram) clearly seen in the 
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Fig. 5. CO(3-2) channel maps of the galaxy NGC 4826. The velocity relative to the systemic velocity of the galaxy (v'h e i=402 km s _1 ) 
is given at the top left of each map. The phase center indicated by the large cross is at a/2000=12 /, 56 m 43.76 I , <5/2000=21 o 40'51.9". 
The galactic nucleus is marked by the small cross. The beam (represented at the bottom left of each map) is 2.6"xl.9" and the rms 
is 75 mJy beam 1 . The only contour line overlaid corresponds to the 4 <x level. 
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Fig. 6. SMA CO(3-2) integrated map of NGC4826 (grey scale given in the color bar in Jy kms -1 beam 1 ) overlaid with the PdBI 
CO(l-O) integrated map (left) with contours at 2, 4, 8, 12, 16, 20, and 24 Jykms -1 beam -1 and the PdBI CO(2-l) integrated map 
(right) with contours at 4.1, 14.3, 24.5, 34.7, 44.9, 55.1, 65.3 and 75.5 Jy kms~' beam -1 . No primary beam correction has been 
applied. All maps have the same resolution (PdBI data have been tapered) and the beam (2.6"xl.9") is shown at the bottom left. 
The white cross labelled C shows the position of the dynamical center (aj2ooo = 12 h 56 m 43*64, <5j2000 = 21°40'59'.'3). The other 
crosses indicate the positions of the spectra plotted in Fig. [8] 
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Fig. 7. Composite color image of NGC4826 made from the intensity maps of the CO(1-0), CO(2-l) and CO(3-2) line emission. 
The line intensities are represented by the red, green and blue color intensities, respectively. The intensity maps were obtained by 
summing all channels where the line has S/N>3 and are corrected for primary beam attenuation. 
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Fig. 8. Individual spectra of the CO(l-O) line (black), the CO(2-l) line (blue) and the CO(3-2) line (red) at 8 different positions 
(see Hig-[6j in the central region of NGC4826 extracted at the same spatial resolution of 2.6"xl.9". The fluxes have been corrected 
for primary beam attenuation and the CO(l-O) flux is scaled up by a factor 3.5 for better readability. 



individual CO(l-O) and CO(2-l) channel maps (INUGA-J) can 
just be recognized in the CO(3-2) channel maps. The integrated 
map (see RGB-image in Fig. [7] and detailed greyscale maps in 
Fig.[6]l, however, shows that the CO(3-2) distribution is similar 
to that of the CO(l-O) and CO(2-l) lines at the same resolution. 
It is dominated by a nearly circular region elongated along the 
galactic major axis with a radius of ~2" (4 pc) and identified 
as a circumnuclear disk (CND) in iNUGA-j . The emission de- 



tected outside the disk seems to coincide with the spiral arms 
described in that same paper. This general agreement between 
the three CO lines suggests that the discrepancies in the channel 
maps are mainly due to the lower S/N of the SMA data in indi- 
vidual channels. This interpretation is further supported by the 
similarities between the individual spectra shown in Fig. [8] At 
the eight positions (marked on the integrated map in Fig.|6]l the 
spectra exhibit similar profiles to within 30%. 
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The relative strengths of the lines appear to follow two 
trends, (i) The closer to the center, the stronger the CO(2-l) 
and CO(3-2) lines relative to CO(l-O). (it) The CO(3-2) line 
is stronger relative to CO(2-l) in the central disk (positions B, 
C and D) and along the major axis (compare positions A and E 
vs F and G), where the spiral arms connect to the disk. 



4. Line ratios 

With the definition adopted in the introduction any line ratio 
would tend to 1 with increasing excitation temperature for ther- 
malized gas in local thermal equilibrium. Line ratio maps were 
computed restricting the spectral window of each pixel to the 
channels where the CO(l-O) flux is >4-cr and the CO(3-2) 
flux is positive. When the spectral window was narrower than 
50km s _1 , i.e. less than 5 channels, the pixel was discarded. 

4.1. NGC4569 

The ratio maps with the corresponding error (standard devia- 
tion) maps are shown in Fig. [9] R21 is distributed in the range 
0.2 to 0.9 with a weighted mean of 0.63. The uncertainties on 
R21 are in the range 0.01 to 0.05. R32 is distributed in the range 

0. 05 to 0.45 with a weighted mean of 0.23. The uncertainties on 
^32 are in the range 0.03 to 0.11. The errors on the weighted 
means based on the data noise only are less than 0.02. While 
calibration uncertainties introduce an overall error of the order 
of 0.05, they should have the same effect on all locations and 
thereby spatial trends should be unaffected. Th e mean values are 
in goo d agreement with the values found by Hafo k & Stutzkil 
(2003) based on single dish observations with a beam of 80" 
HPBW, i.e. fl 2 i =0- 64+0.08 and i ?32=0. 23±0.04; and with the 
R32 value found by IWilson et alj d2009l) within the inner 22", 

1. e., 0.34+0.11 (note however that the latter is not a fully inde- 
pendent check as we used the same CO(3-2) single dish data to 
complement our SMA data). 

The R21 and R22 maps are different. While ^21 is smooth, 
more or less symmetric with respect to the center, and continu- 
ously increasing toward the center with a maximum at ~\" to 
the west of the galactic center, 7?32 is more clumpy, exhibiting 
several peaks, and i s more asymmetric. The large-scale gradient 
in /?32 reported by Wilson et al. (2009) does not appear in the 
inner 10". 

R12 is expected to probe more extreme conditions such as 
those occuring in star forming or shock regions. This is con- 
firmed by the comparison to the Paa emission (FigurefTUb. 
Although the detailed structure of ^32 is not reflected in Paa, 
most of the pixels where 7?32>0.2 coincide with Paa emission de- 
tected in the NICMOS image. Rz\, instead, reaches high values 
outside (to the north) as well as inside the Paa emission region, 
although it peaks at the galactic center as the Paa. The fact that 
R12 tracks massive star formation in NGC 4569 better t han /?2i 
is consistent with the conclusion o f IWilson et alJ {2009), based 
on maps over larger areas, that CO(3-2) is a robust tracer of star 
formation. 

The fact that R^ correlates with the individual CO line inten- 
sities (the three peaks in CO are also peaks in R22) also suggests 
that it is more sensitive to density. In contrast, the fact that /?2i 
does not show the same structures as the individual line maps 
and that it increases steadily towards the center suggests that it 
is more sensitive to temperature. This is expected for thermal- 
ized gas with moderate opacity in the CO(3-2) line. 



The left panel of Fig. Q~T] shows the distribution of the ratios 
in the R12 vs R2\ space; in the following this will be referred to 
as the /?32-vs-/?2i diagram. This diagram was obtained by map- 
ping each pixel (0.2" xO. 2" in size) of the ratio maps to the grid 
and convolving with a gaussian kernel of 0.02 FWHM, corre- 
sponding to a typical uncertainty in the line ratio. Based on the 
line ratio error maps shown in Fig. [9] it should be kept in mind 
that any feature in this diagram with a size lower than ~0.05 can- 
not be considered as real. In addition, although the ratio maps 
are considered as continuous distributions (we used small pixels 
with respect to the beam), it should be kept in mind that the reso- 
lution is finite in the ratio maps. This diagram is used to show in 
a synoptic way which regions of the ratio-ratio space are popu- 
lated. As expected from the differences between the R32 and /?2i 
maps noted above, the distribution is dispersed in the diagram. 
The ratios populate a region going from (0.4, 0. 15) to (0.8, 0.25) 
with an average thickness of ~0.15 in an orthogonal direction. 

4.2. NGC 4826 

The ratio maps with the corresponding error (standard deviation) 
maps are shown in Fig. [12] 7?2i is distributed in the range 0.2 
to 0.95 with a weighted mean of 0.67. R^2 is distributed in the 
range 0.05 to 0.60 with a weighted mean of 0.38. The errors 
on the ratios and the weighted means are similar to those in the 
case of NGC 45 69. Th e mean values are lower than the ratios 
found bv llsraell (120091) based on single dish observations in the 
inner 21", i.e. fl 2 i=0.98+0.17 and fl 32 =0.59+0.09 and the in- 
ner 12", i.e. /?2i = 1.06+0.22 and /? 32 =0.79±0.18. They are, how- 
ever, within (or close to) the 2<x interval of the single dish mea- 
surements, does not appear to be more clumpy than R21, in 
contrast to NGC 4569, and both maps look relatively symmetric 
with respect to the dynamical center, with both ratios increasing 
steadily toward the center. There are however differences. While 
the 7?2i spatial distribution appears to be close to circularly sym- 
metric, the /?32 distribution is more elongated along the galactic 
major axis. Furthermore, R21 seems to be more prominent to the 
north of the dynamical center (although its peak coincides with 
the dynamical center), whereas R22 forms a ridge along the ma- 
jor axis south of the dynamical center with its peak ~1" to the 
southwest. 

As in NGC 4569, high R i2 values (> 0.3) occur within the 
Paa emission region as detected by NICMOS, while this is not 
necessarily true for R21 (see Fig. Qj] right panel). Indeed, al- 
though R21 increases toward the center as R32, it reaches high 
values (>0.7) north and south of the Paa emission region, which 
is elongated along the galactic major axis. R21 shows some cor- 
relations with the 2 cm radio continuum (see Fig.[13]left panel). 
In particular, the 2 cm peaks 1" and 4" north of the dynamical 
center coincide with high R21 values. At the VLA's resolution, 
the 2 cm continuum is dominated by emission from HII regions 
and SNRs around which the molecular phase is expected to be 
strongly heated. If R21 is more sensitive to temperature relative to 
R32, this may explain the better correlation of 7?2i with the 2 cm 
continuum than R^- We emphasize, however, that this trend is 
weak relative to the radial increase toward the center of R 21, R32 
and P aa emission (as for R21/R32 , previously noted by Israel 
2009). In fact both line ratios are sensitive to temperature and 
density and their interpretation in terms of physical conditions 
requires radiative transfer modeling; this is the topic of the next 
section. 

The middle panel of Fig. QT| shows the /?32-vs-/?2i diagram 
for NGC 4826. The emission is much more widely spread in the 
^32-vs-/?2i diagram than for NGC 4569. The top part of the dia- 
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Fig. 9. Line ratio maps of NGC4569 and associated uncertainties. Upper Left: R21 with black contours ranging from 0.3 to 0.7 
in steps of 0.1 and white contours at 0.8 and 0.9. The beam (2.3"x2.1") is represented as a filled ellipse at the bottom left corner. 
Upper Right: R32 with black contours ranging from 0.1 to 0.25 in steps of 0.05, and white contours at 0.3, 0.35 and 0.4. Lower 
Left: R21 error map with black contours at 0.02, 0.03, and white contours at 0.04, 0.05. Lower Right: R32 error map with black 
contours 0.05, 0.07, and white contours at 0.09, 0. 1 1, 0. 13. 



gram with R32 > 0.3 being more populated, in particular the top 
right part with R 2 i > 0.6 and R32 > 0.3. 

5. Modeling the physical conditions of the 
molecular gas 

5. 1 . Pseudo-LTE model 

The line intensity of a single component medium in local ther- 
mal equilibrium (LTE) is entirely determined by the tempera- 
ture of the medium, T, and by the ratio of the column density 



to the local line width, Nco/dV when the line is not optically 
thick. In reality, the ISM is multi-component and the intensity of 
each line results from the sum over the beam of all the individ- 
ual components' line intensities. LTE may still be a reasonable 
approximation for the local emission of individual components 
or clumps, but the average opacity in a line toward the observer 
also depends on more global properties of the ISM (e.g. turbu- 
lence and complex structures). A way to improve the model is 
to decouple the average opacity from the average local physical 
conditions and consider a mixture of components. For the sake of 
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Fig. 10. Line ratio maps with the same contours as in Fig.|9](only above 0.2 for ^32) overlaid on the Paa map obtained from NICMOS 
data (Programme ID: 9360, PI: Kennicutt). The grey scale is logarithmic; dark color corresponds to high flux. 



simplicity we assume that CO(l-O) and CO(2-l) emissions are 
optically thick, and we consider optically thick (t>>1) and op- 
tically thin (t«1) components for the CO(3-2) emission. This 
choice is justified by the variations of the line opacities with the 
physical conditions as shown in Fig. [14] and based on large ve- 
locity gradient model (LVG model described in Section [5^2b . For 
gas with a temperature close to or below ~ 20 K and a density 
above ~ 10 3 cirT 3 , the CO(l-O) and (2-1) emission is mainly 
optically thick whereas the CO(3-2) emission optical depth can 
vary from thick to thin. We denote fa the fraction of optically 
thin emission for CO(3-2). Treating fa as an independent pa- 
rameter allows us to avoid complex modeling of the ISM, but it 
should be kept in mind that this parameter is also related to the 
physics of the ISM. In other words, the interpretation of such 
a model requires further physical modeling (at least qualitative) 
of the opacities. We assume the line intensities are given by the 
LTE equations (the optically thick component is independent of 
Nco/dV). We refer to this modified LTE model as a pseudo-LTE 
model. Although this model relies on simple assumptions, we 
think it can be used to interpret the overall distribution of the 
gas in the R32-VS-R21 diagram in terms of variations of its physi- 
cal properties. In practice, we reproduce the overall distribution 
in the R32-VS-R21 plot with families of pseudo-LTE models cor- 
responding to different combinations of physical parameters. In 
other words, this model is not used to fit each (R21, ^32) pair 
independently (this would not be possible given the number of 
unknowns); it is used to relate the different parts of the ^32-vs- 
R21 diagram to the physical conditions and reproduce the overall 
ratio distribution. 



5.1.1. NGC4569 

As shown in the right panel of Fig.Qj] a single-component 
LTE model (curves 1 and 2) cannot reproduce the line ra- 
tios observed in NGC4569 (cyan contour). Even with high 
column densities per unit velocity (curve 1 corresponds to 
Nco/dV = lO^cirT^kms- 1 )- 1 and curve 2 to N C0 /dV = 
12 x 10 15 crrT 2 (kms^'f 1 ) at the low-temperature (lower left) 
end 7?32 is always lower than observed. 

Curves 3 and 4 correspond to pseudo-LTE models with 
/32=0.5 and fa=0.2 respectively. The column density to line 
width ratio is Nco/dV = 10 15 cm~ 2 (kms -1 ) -1 and the temper- 
ature goes from 4 to 20 K along the curves from bottom left 
to top right. Neither curve can fit the overall ratio distribution, 
although each one intersects with the distribution. The ratios 
corresponding to the bottom left of the diagram correspond to 
gas at ~5 K with the three lines being mainly optically thick 
(curve 4) and the ridge at the top right of the distribution to 
gas at -15 K with 50% of the CO(3-2) line emitted by optically 
thin gas (curve 3). Thus, the overall distribution of the ratios 
in the /?32-vs-7?2i diagram noted in the previous section could 
correspond to gas spanning the whole range of physical proper- 
ties between these two limits, i.e. (T-5, /32=0.2) and (T=15 K, 
/32=0.5). For the gas to remain within the observed ratio dis- 
tribution the CO(3-2) thin fraction, fa, must increase with the 
temperature. To qualitatively interpret the variations of the op- 
tically thin fractions we consider the opacities given by a large 
velocity gradient model for a range of temperatures and densi- 
ties and displayed in Fig. [14] Qualitatively, gas that is optically 
thick in the three lines is denser (n > 10 4 otT 3 ) than gas that 
is optically thin in the CO(3-2) line but still thick in the other 
two lines (10 3 < n < 10 4 cm 3 ) for T<20K. Hence, increas- 
ing fa means increasing the fraction of more diffuse gas, and 
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Fig. 11. (Left) 7?32-vs-/?2i diagram for the galaxy NGC4569. The x-axis represents R21 and the y-axis ^32. The intensity scale gives 
the sky area in the map per unit area in the ratio space, it is expressed in beams per 0.01 x 0.01 ratio area. This diagram has been 
produced by convolving the actual distribution of pixels in the ratio-ratio plane by a gaussian kernel with a FWHM of 0.02. Note 
that the pixels in the ratio maps used to compute the distribution are much smaller than the beam (see Section[4]i. The contour levels 
are 0.3, 0.6, 0.9, 1.2 and 1.5 beams/0.01 2 . (Middle) ^ 32 -vs-/?2i diagram for the galaxy NGC4826; the contour levels are 0.3, 0.5, 
0.7, 0.9, 1.1 beams/0.01 2 . (Right) LTE and pseudo-LTE models represented as black curves in the ^32-vs-/?2i diagram. Curves 
1 and 2 correspond to pure LTE models with Nco/dV = 10 15 cirr 2 pc (kms -1 ) -1 and Nco/dV — 12 x 10 15 cm~ 2 pc (kms -1 ) -1 , 
respectively. The temperature increases from 4 K at the bottom left to 8 K at the top right. Curves 3 and 4 correspond to pseudo-LTE 
models with / 3 2 = 0.5 and fa = 0.2 respectively (see Section lBTTI for the definition of fa)- The column density to line width ratio 
is Nco/dV = 10 15 cirT 2 pc (km s -1 ) -1 and the temperature increases from 4 to 20 K along the curves from bottom left to top right. 
The cyan and magenta curves show the lowest contours for NGC 4569 and NGC 4826, respectively. 



from left to right the ratio distribution (i.e. the distribution in the 
7?32-vs-/?2i diagram) corresponds to increasing temperature and 
decreasing average density. 

5.1.2. NGC4826 

As shown in the right panel of Fig.Q~T| The highest ratios ob- 
served in NGC 4826 (magenta contour) follow more or less the 
pseudo-LTE model represented by the curve 4 in the 7?32-vs-/?2i 
diagram. As noted in Section l4~2l the line ratios increase toward 
the galactic center (FigfT2l. although their maxima do not ex- 
actly coincide with it. Hence according to this simple model the 
temperature should increase toward the center and a significant 
fraction (~20%) of the CO(3-2) emission should be optically 
thin. 

5.2. LVG model 
5.2.1. The procedure 

To interpret the observations of the three CO lines in terms of 
molecular gas physical conditions we have modeled the line 
ratios based on the large velocity gradient (LVG) a ssumption 
(IScoville & Solomonl[T97llGoldreich & Kwanlll974l) using the 
LVG code in the MIRIAD package. As opposed to the LTE as- 
sumption, the LVG assumption allows us to deal with transitions 
that are not necessarily thermalized. On the other hand, with 
the observations available we can only consider a single compo- 
nent. Although the single component assumption is at odds with 
the known complexity of the real ISM (even at the 40 pc lin- 
ear resolution reached in NGC 4826), it was successfully used 
in many extragalactic studies to assess the physical properties 
of the molecular gas. For a single component the free parame- 
ters are the molecular hydrogen number density, n(H2), the ki- 
netic temperature, T^, and the ratio of the CO molecule frac- 



tional abundance with the velocity gradient, Xco/(dV/dR). With 
only two observables (the line ratios) one of the three parameters 
needs to be fixed. 

We have implemented a maximum likelihood fitting proce- 
dure to fit two parameters assuming normal probability distribu- 
tions for the line ratios. The normalized likelihood of the best fit 
is used as a goodness-of-fit metric, denoted Al. It is equal to 1 
when the observed line ratios are exactly reproduced by the best 
fit model and it tends to with increasing discrepancy. 

The value of Xco /(dV/dR) was constrained as follows. We 
have generated three grids of models for Xco/(dV/dR)=10~ 6 , 
10~ 7 and lO^pc^kms -1 ) -1 . Each grid is 50x50 in size corre- 
sponding to 50 values of 7\ in the range 5 to 200 K and 50 val- 
ues of n(H2) in the range 10 2 to 10 5 cm 3 . The intervals in 7\ and 
n(H2) are logarithmic. The grids are displayed in Fig. [15] To es- 
timate the most appropriate value of Xco /(dV/dR) we computed 
Al for each point of a R32-VS-R21 grid assuming an uncertainty 
of 0.02 for R 2 i and 0.04 for R 32 - The distribution of Al over the 
^32-vs-/?2i grid thus obtained was then compared to the ^32-vs- 
/?2i diagram of the galaxies. The three Al grids are plotted in 
Fig-E] For both galaxies the Xco/(dV/dR) value maximizing 
the goodness-o f-fit over the ma p is 10 _7 pc(kms _I ) 1 similar to 
that found bv llono et all (120071) for NGC 6240 following a simi- 
lar method; we therefore adopted this value for all the models. 

As can be seen in Fig. [15] the two line ratios follow a sim- 
ilar distribution over the density-temperature grid, i.e. the iso- 
contours of both line ratios are close to parallel over most of the 
grid. Hence, in many cases, when the uncertainties are taken into 
account the solution is degenerate. The confidence intervals typ- 
ically have arc-shapes similar to the isocontours (see Fig.[T7land 
[191. However, a confidence interval does not cover the complete 
range of densities and temperatures, and the average values over 
this interval can be used to characterize qualitatively the trend 
of the model, i.e. whether it tends to favour high or low densi- 
ties and temperatures. Thus, instead of considering the best fit 
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Fig. 12. Line ratio maps of NGC 4826 and associated uncertainties. The crosses are the same as in Fig. [6] Upper left: R21 with black 
contours at 0.4, 0.5, 0.6, 0.7, 0.8 and white contours at 0.77 and 0.9. The beam (2.6"xl.9") is represented as a filled ellipse at the 
bottom left corner. Upper right: R32 with black contours at 0.1, 0.2, 0.3 and white contours at 0.4 and 0.5. Lower left: R21 error 
map with black contours at 0.02, 0.03 and white contours at 0.04 and 0.05. Lower right: R32 error map with black contours at 0.05, 
0.07 and white contours at 0.09 and 0.1 1 



solutions, we consider the average values of the temperature and 
density over all possible solutions, i.e. the average values over 
the 90% confidence interval. As the 90% confidence interval can 
be unbounded toward high temperatures and/or high densities 
we fix the maximum possible temperature to 70 K and the max- 
imum possible density to 10 s crrT 3 . Fixing these upper limits 
affects the averages but not their qualitative variations with the 
ratio values. 



5.2.2. NGC 4569 

The LVG fits to the eight individual positions marked on the map 
of Fig. |2] are displayed in Fig. [17] This figure shows that good 
LVG solutions exist for all positions but that the uncertainties are 
large (see the confidence intervals in the insets). This figure also 
shows how using the average values instead of the best fit gives 



more stable results. For example, the ratio values at the northern 
and southern peaks (B and G) are similar and should therefore 
give similar results. The best fit solutions are very different, the 
average values are more stable and therefore better adapted to a 
qualitative description. 

Density and temperature maps obtained by applying the fit- 
ting procedure to all pixels are presented in Fig. [18] The pres- 
sure, obtained by averaging the product, T x n over the confi- 
dence interval, is displayed at the bottom left and the goodness- 
of-fit map is shown at the bottom right. Gas with density above 
5.5 x 10 3 cm -3 (the lowest white contour) and temperature below 
15 K is distributed around the center in an elongated ring extend- 
ing from the CO(1-0) northern peak (position B) to the CO(1-0) 
southern peak (posi tion G). Th is is consistent with the kinematic 
model presented in NU GA- VTJl where the molecular gas was dis- 
tributed in an elongated ring extending from the northern peak 
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Fig. 13. Map s of N GC4826. Left: /?2i with the same contours as in Fig. [12] overlaid on the 2cm r adio conti nuum map from 
iTurner & Hoi dl994l) . Right: R 32 with the same contours as in Fig. [12] overlaid on the Paa map from iNUGA-fl Dark color cor- 
responds to high flux. 
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Fig. 14. Variations of the CO line opacities with respect to temperature and density based on LVG modeling. From left to right 
panel: CO(l-O), (2-1) and (3-2) opacities. Used for the qualitative interpretation of the pseudo-LTE model. 



to the southern peak and identified as the ILR. Within this ring 
there are two clumps of high density cold gas, one at the south- 
ern CO(l-O) peak (position G) and the other 3" to the west of 
the nucleus (position F). 

In the central region delineated by the ring there is also a den- 
sity peak ~1.5" north of the center (near position D) where the 
temperature is higher than in the ring. On average, in the central 
region the density is higher than outside the ring. The tempera- 
ture is higher in the central region than in the ring itself but is 
not higher than in the region outside the ring. This reversal is not 
expected, as the central starburst should heat the gas, and this 
may reflect the limit of the method used here, i.e. averaging the 
physical parameters over the confidence interval can only give a 
trend not the absolute values. The warmest regions are to the east 
of the northern peak and to the south-west of the southern peak. 
They may correspond to diffuse hot gas near the points w here 
the orbit crowding causes shocks (see orbits in lNUGA-VIlh . We 



note however that south of the southern peak the goodness-of-fit 
is lower (Ai < 0.4). There are three pressure peaks coinciding 
with the three density peaks. 



5.2.3. NGC4826 

The LVG fits to the eight positions marked in the map of Fig. [6] 
are displayed in Fig. [19] This figure shows that good LVG solu- 
tions exist for all positions except the easternmost and western- 
most ones (positions A and E) where R32 is too high with respect 
to #21. 

Density and temperature maps obtained by applying the fit- 
ting procedure to all pixels are presented in Fig. [20] The pressure 
is displayed at the bottom left, and the goodness-of-fit map is 
shown at the bottom right. The densest (with density > 10 4 cm -3 ) 
and coldest (T ~ 5K) gas is distributed in a semicircle ~60- 
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Fig. 15. Grids of LVG models for a range of molecular hydrogen densities and kinetic temperatures. Each grid corresponds to 
a different value of X C o/(dV/dR), from left to right: 10" 6 , 1(T 7 and 10 _8 pc(kms _1 ) _1 . The CO(l-O) radiative temperature is 
represented in color scale and is overlaid with R21 black contours (0.3 to 1.3 in steps of 0.2) and blue contours (0.1 to 0.6 in 
steps of 0.1). The line ratios increase from bottom left to top right. 
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Fig. 16. Likelihood of all possible (RiuRn) pairs on a /?32-vs-/?2i grid for different values of Xco/(dV/dR). The contours of the R32- 
vs-/?2i diagrams of NGC 4569 (blue) and NGC 4826 (magenta) are overlaid. The contours are smoothed and their levels correspond 
to 3 and 9 pixels per grid bin. (see Fig.fTTTi. 



80 pc in radius whose center coincides with the dynamical cen- 
ter. Inside the region delineated by the arc the temperature in- 
creases from the arc to the north of the dynamical center. 

However, as can be seen from the goodness-of-fit map in the 
lower right panel of Fig.|20]and as was already noted from the fit 
to the individual positions A and E (Fig.[T9l>, the LVG solutions 
outside the arc along the major axis, i.e. the solutions giving the 
most diffuse and warm conditions, do not fit the observed ratios 
well. 



6. Discussion 

6.1. Comparison of the two galaxies 

NGC 4569 and NGC 4826 both harbor transition type nuclei. It 
is therefore interesting to compare the ISM close to the nucleus 
and look for similarities. At first sight, the molecular gas dis- 
tributions and properties appear to be very different in the two 
galaxies. In NGC 4569, the gas is concentrated in an ILR ring 
and a clump close to the center, and the CO is not very excited. 
In NGC 4826, the gas is concentrated in a CND and the gas is 
more excited. In addition, the strong correlation between the two 
line ratios near the center in NGC 4826 is not seen in NGC 4569. 




Fig. 21. Pressure contours of NGC 4826 (same contours as in 
Fig. [20] lower left panel) overlaid on the Paa map. The ellipse 
corresponds to a ring of 90 pc radius inclined by 60deg. The 
straight line indicates the major axis of the galaxy. 
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Fig. 17. LVG best fits to the line ratios measured in the central region of NGC4569 at the eight positions marked in Fig. [2] The 
measured line ratios are represented with the 1-cr error bars and the best fit model ratios are linked with a straight line. The inset 
in each plot shows the 90% confidence interval (light grey shaded area) and the 50% confidence interval (dark grey) in the density- 
temperature plane assuming a normal probability distribution for each ratio; the best fit solution is marked with a cross. The best 
fit temperature and density as well as the averages over the 90% confidence interval are given at the top right of each plot. The 
normalized likelihood of the best fit, P/P max is also given (1 means the measured values are exactly reproduced by the model). 
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Fig. 18. LVG result maps for NGC4569. Upper left: H2 number density averaged over the 90% confidence interval of each LVG 
model. The black contours are at 5 x 10 2 , 10 3 , 3 x 10 3 cm 3 and the white contours at 5.5 x 10 3 , 10 4 , and 5 x 10 4 cm 3 . Upper right: 
kinetic temperature averaged over the 90% confidence interval of each model. The black contours are at 5, 10, 15, 20, 25 K and the 
white contours at 30, 35, 40 K. Lower left: pressure averaged over the 90% confidence interval of each model. The black contours 
are at 10 , 10 44 , and 10 47 KcirT 3 and the white contours at 10 5 , 10 5 3 and 10 5 6 Kcm ~ 3 . Lower right: the normalized likelihood 
(or merit), P/P max , of the best fit with black contours at 0.35 and 0.55 and white contours at 0.75 and 0.95. 



However, the low excitation gas in NGC 4826 (with /?32<0.3) 
follows a ratio distribution similar to that of NGC 4569 (see 
the /?32-vs-7?2i diagrams Fig.fTTTi. and the difference in linear 
resolution must be taken into account. It is 4 times higher for 
NGC 4826 (40pc) than for NGC 4569 (160pc). If NGC 4826 
were observed with the same linear resolution as NGC 4569 the 



CND would not be resolved and the ratio distribution in the up- 
per right of the R32-VS-R21 diagram would be averaged into a 
small clump in the diagram. Conversely, observing NGC 4569 
with the same linear resolution as NGC 4826 could well reveal 
more excited gas. 



18 



F. Boone- et al.: Mapping the physical conditions of the molecular gas in the central regions of two active galaxies 



FAR EAST 



n.i 



T. = f E u-=1.3n+Ci 



1 



C) 



-■■i-Vi-D/v-y.l-w- i.^:»-aj/r:t>(i-(ij 
DYNAMICAL CENTER 





T,= HE !.=:■!. Se +01 riii-* 










T. = 11 E n =: : ..-": -KM m": 










it? 


















+ 








1 















E) 



FAR WEST 



coru-al/raci-a) 



D.f! 



!).3 





. T,= - 

r /!■__■ 

T. = 


E i_=J.£i!!-KS-i fui"? 

' (i.tlfl | 

" E ii. =t).Gr1-Ul m" : "j 






V 














1 







G) 



"Ois-lj/t'oii-'J) 

SOUTH 



0.7 r 



T,= 9 E iv=4.Hc4ai cm*" 

i : iw= |a« 

T, ■ 5 E .li cm" 



u.l 

f'ote-ij/cau-n) 



B) 



EAST 



1 1V1 : __= fcfiffl 



a 1 



D) 



rO-i'j-lVVua-U) 
WEST 



[.■&( :s— y i /■"!>( 1— c i 



U.7 



D.O 



T,= iLi E Ii-=U..at+S.l otcf? 
1" 1„F 0.9*J 

T. = U E li al5.6.-h04 0& 



4 



F) 



NORTH 



(.'(?( :i-i:)/am^(l) 



0.3 



I P /!'__■ i.uy 



e ui.Vr+ea fin 



T. = 1?> E n - = l.lc-+-i>:i in" 



IG 1 



.£ 1 



+ 



3 1 5 



D.2 



CQta-jj-j/cufl-m 



H) 



FAR SOUTH 



d.1 



0.!J 



T,= Bii K il=.J.Hc+SJ 'ti"' 
P/l l -aF ^.36 
T, - 3-1 E n -il.l: +ij:> ti^ 



\ 



n.i 

CO^-l)/t'Dl|i- D J 



? :> 1 3 



Fig. 19. LVG best fits to the line ratios measured in the central region of NGC4826 at the eight positions marked in Fig. [6] See 
caption of Fig.[17]for more details. 
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Fig. 20. LVG result maps for NGC4826. Upper left: H2 number density averaged over the 90% confidence interval of each LVG 
model. The black contours are at 10 3 and 5 x 10 3 cirT 3 and the white contours at 10 4 and 5 x 10 4 cirT 3 . Upper right: kinetic 
temperature averaged over the 90% confidence interval of each model. The black contours are at 5, 10, 20 K and the white contours 
at 40 and 60 K. Lower left: pressure averaged over the 90% confidence interval of each model. The black contours are at 10 47 , 
10 49 , 10 51 and 10 3 KcirT 3 and the white contours at 10 5 5 and 10 5 7 KcirT 3 . Lower right: goodness-of-fit with the same contours 
as in Fig. [18] 



Taking into account this observational bias, the most robust 
similarity between the two galaxies is the relatively low excita- 
tion of the CO compared to other galax y centers. At the hundred 
parsec scale, R?\ and Ry can be >0.8 feckart et al.lfl990l 1 199 it 



parsec scaie, K21 anu K32 can De >u.o (ncK.ari ei ai. ryyu, ivyi; 
Braine & Combeslll992t IWild et aljfl992l) and >1 dDumke et al.1 



2001), respectively, in galaxy centers with enhanced star for- 
mation, and the ratios lie typ ically in the range 0.7—1.1 and 
0.5-1 around luminous A GNs (Papadopoulos & Seaquist 1999; 
Papadopo ulos et ail 120081: Ilsraelll2009l) . At the same scale ^21 
even reaches a value of ~2 in the p rototypical Seyfert galaxy 
NGC 1068 dSchinnerer etal.ll200qbl) as w ell as in the Seyfert 



galaxy NGC 1097 dHsieh et al.1 12008). R32 also reaches 



value of 2 at the hundr ed parsec scale in the nucleus of M51 
(Matsushita et ai] |2004l) . Finally, R32 is also lower than in the 
inner 200 pc of the Milky Way (a.k.a. the circumnuclear zone, 



CNZ), where R y has an average value of ~0.8 (see Fig. 3 of 
lOka et all2007l) . Thus, at the scale of previous single dish stud- 
ies and at the scale of the CNZ, the molecular gas in the centers 
of NGC 4569 and NGC 4826 appears to be in a rather low excita- 
ti on state. R21 and R32 are al so lower than the model predictions 
of IWada & Tomisakal d2005l) for AGN tori. 



6.2. Can opacity effects explain the asymmetries in 
NGC 4826? 

The LVG modeling indicates the presence of a semicircular arc 
of dense gas close to the nucleus of NGC 4826. Its center of cur- 
vature agrees with the dynamical center and it lies at the south- 
west of the nucleus. In fact, this dense arc seems to reflect asym- 
metries in the CO and R32 maps (Figs. l6l and [T2l. Indeed, close 
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Fig. 22. Sketches illustrating the two opacity effects considered in Section Kx2l In both panels the color scale is related to the 
temperature in K as shown in the colorbars and the star symbol represents the AGN. The observer is assumed to be looking from 
the top right, so that the right-hand side of the panels corresponds to the near side and the left-hand side to the far side of the galaxy. 
The top panel illustrates how the densest and coldest gas distributed in a ring at 70 pc from the center could partly obscure emission 
from the dense torus on the near side and from the galactic disk on the far side. The white contours correspond to density levels. The 
dashed white lines delineate the regions that would be obscured by the densest gas. The lower panel illustrates how clumps heated 
by the nucleus could have inhomogeneous properties with the side exposed to the nucleus being warmer, less dense and less opaque 
to CO emission than the other side. 



to the nucleus the emission of the three CO lines is slightly more 
prominent south of the nucleus and, as was noted in Section l4~2l 
the /?32 distribution is clearly more prominent to the south. 

Interestingly, as can be seen in Fig.[2T] the dense arc seems 
to coincide with regions of the CND where the Paa emission 
is weaker. Another interesting point is that the galaxy major axis 
seems to be a privileged axis for this anticorrelation, i.e., the Paa 
emission is more prominent on the near side (to the north east) of 
the galactic disk, whereas the semicircular dense arc is on the far 
side; the minor axis seems to be a symmetry axis for bo th the Paa 
emission and the dense gas distribution. As noted in INUCiA-1 
the dust lane extinction is usually stronger on the near side (an 
effect that turns out to be spectacular in this galaxy outside the 
CND), and the Paa near side/far side asymmetry is opposite to 
the one expected for star formation in a density wave. A possi- 
ble explanation for this Paa asymmetry is that extinction is not 



confined to dust lanes in a density wave but rather comes from 
extraplanar dust around the nucleus (see e .g. iMason et"ail 2006: 
lRhee &Larkinll2006t lOuillen et al.ll2008l) . In addition, part of 
the Paa emission close to the nucleus may come from material 
above or at the surface of molecular clouds that has been ionized 
by the nuclear activity rather than by young stars distributed in a 
spiral. 

The near side/far side asymmetries and the Paa/dense gas 
anticorrelation suggest that the CO emission as well could be af- 
fected by opacity (for a similar conclusion b ased on int ensity ra- 
tios averaged over a larger beam, see also llsrae 3120091) . Opacity 
effects are generally not considered for extragalactic CO obser- 
vations. Although the low-/ CO lines are mainly optically thick, 
the dumpiness of the ISM and the very low probability to have 
two clumps on the same line of sight with the same line of sight 
velocity imply that molecular clouds do not shadow each other, 
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CO(3-2) and CO (1-0) contours 
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Fig. 23. The upper panels show the CO integrated maps produced by the flared disk toy model described in Appendix lA.il The 
CO(3-2) emission is shown in greyscale, and the CO(l-O) (left panel) and CO(2-l) emission (right panel) are contours. The lower 
left and right panels show the resulting /?2i and R32 maps, respectively. 



and t herefore that CO line s are good column density probes (see 
e.g. iDickman et al.llT986l) . CO morphologies observed in galax- 
ies are thus assumed to reflect the actual gas distribution. We see 
however two possible reasons for opacity to play a role in the 
CO morphologies observed close to galactic centers. First, the 
column density is larger and the probability of having several 
clumps on the same line of sight with the same line of sight ve- 
locity may not be negligible anymore. Second, the activity of the 
nucleus is expected to affect the surrounding ISM and as a result 
the individual clumps may have strongly inhomogeneous prop- 
erties fully determined by their distance and orientation relative 
to the nucleus; the clumps on the far side of the galaxy could 
thus show very different properties from those of the near side. 
Both effects, which we refer to as the large scale opacity effect 
and the clump scale opacity effect, respectively, are sketched in 
Fig. [22] We investigate them in more detail in the following. 



6.2.1 . Large scale opacity effect 

The molecular gas disks in active nuclei are thought to flare 
out to radii of order lOpc, forming the putative obsc uring 
torus required by the unified AGN model (see e.g. lElitzurl 
2008, and references therein). How this structure connects to 
the outer disk and whether it is embedded in a larger toroidal 
structure is still uncertain . The hydrodynamical simulations by 
IWada & Tomisakal (2005) suggest that such a large molecular 
structure could extend out to several tens of parsecs and this 
seems to be confirmed by observatio ns of the near-IR m olecu- 
lar hydrogen lines in nearby AGNs by Hick s et a"D (I2009I) . 

If the densest and coldest gas is distributed in a ring inside 
such a large toroidal structure and is able to partly absorb CO 
emission along the line of sight, then the disk beyond the ring 
would be partly absorbed on the far side and the warmer inner 
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Fig. 24. The upper panels show the CO integrated maps produced by the clumpy flared disk toy model described in Appendix lA.2l 
The CO(3-2) emission is shown in greyscale, and the CO(l-O) (left panel) and CO(2-l) emission (right panel) are contours. The 
lower left and right panels show the resulting Rj\ and 7?32 maps, respectively. 



regions of the torus would be partly absorbed on the near side 
(see sketch in the top panel of Fig.1221. In other words, the line 
of sight toward the torus on the near side would be more contam- 
inated by emission from the disk and would not be seeing all the 
emission from the inner regions of the torus. The torus would be 
more clearly visible on the far side. 

In order to check whether this effect could be compatible 
with the observations we consider a toy model of a flared disk 
(see details in Appendix lA.ll i. The resulting integrated CO maps 
and line ratio maps are presented in Fig. [23] Although this simple 
model cannot reproduce exactly the integrated maps shown in 
Fig-El it can partly reproduce the morphology inside the CND. 
In particular, the emission does not peak at the nucleus but to 
the east and west of the nucleus, and there is more emission to 
the south (the far side) than to the north of the nucleus. Also, the 



general morphology does not change significantly in the three 
CO lines. The main differences between the /?2i and maps 
noted in Section l4~2l are also reproduced by the model. The R32 
distribution is more extended along the major axis and is more 
prominent to the south of the nucleus (on the far side) than the 
R21 distribution. 

6.2.2. Clump scale opacity effect 

Alternatively, the asymmetry could result from the inhomogene- 
ity of the clumps forming the toroidal structure. The side of the 
clumps exposed to the nucleus would be warmer and more dif- 
fuse than the other side. The dense side would be more opaque 
and it would hide the warmer side to the observer depending 
on the angle between the line of sight and the line subtended 
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by the clump and the nucleus (see sketch in the lower panel 
of Fig.l22l). A similar "external heati ng" scenario was already 
discussed bvlBaker & Scovilld (1 19981) and lBakerl (120001) and by 
iMeier & Turner! (1200 ll) to interpret asymmetries of molecular 
line emission in NGC 1068 and IC342, respec tively. A simi- 
lar arg ument was also invoked and formalized bv lNenkova et al.l 
(2008b a) to interpret the infrared observations of AGNs. To test 
whether such an effect could be consistent with the observations 
we have distributed clumps in the same toroidal structure as the 
one used in the model of the previous subsection. The resulting 
CO maps as well as the ratio maps are presented in Fig. [24] They 
show that this simple model can also reproduce the slight CO 
emission excess on the far side and the differences between the 
two ratio maps noted in Section l4~2l 

6.2.3. A common property of active galaxies? 

We note that in NGC 4569 as well, the CO emission and the 
i?32 distributions close to the nucleus are more prominent on 
the far side (east side). In fact, although this remained unno- 
ticed for most of them, several nearby active galaxies show a 
slight near side/far side asymmetry close to their nuclei with 
CO emission excess on the fa r side as well; this is for exam- 
ple the case for Centaurus A (Espad a et al.ll2009l) . NGC 4945 
(IChou et all 120071). NGC 3227 dSchinnerer et alf l^OOa). and 
NGC 1068 Jschinnerer et al.ll2000bl) . A detailed statistical study 
of CO asymmetries in nearby AGNs is beyond the scope of this 
paper, but these few cases suggest that near side/far side asym- 
metries could be common in AGNs. 

While both torus models described above can reproduce the 
asymmetry observed, the first model requires a high inclina- 
tion, a very specific geometry, and a high volume filling factor 
of the dense gas. The second model is more robust to various 
inclinations (although it also decreases with decreasing incli- 
nation) and configurations; it would therefore provide the best 
explanation for a putative generalized near side/far side asym- 
metry in AGN s . An alternative opacity effect was proposed by 
iGalliano etal] d2003l) to explain the asymmetry in NGC 1068. 
They invoked CO line absorption by a thick disk, and in their 
model the asymmetry resulted from a tilt between the molecu- 
lar disk and the Compton thick absorber. This model however 
does not explain why the asymmetry should coincide with a far 
side/near side asymmetry. 



7. Conclusion 

We have observed the emission from the first three 12 CO /-lines 
in the nuclei of two nearby galaxies, NGC 4569 and NGC 4826, 
with the Plateau de Bure interferometer and the SubMillimeter 
Array. Combining these data allows us to compare the emission 
in the three lines and to map the line ratios R21 and R32 at a res- 
olution of ~2", i.e. a linear resolution of 160 pc for NGC 4569 
and 40 pc for NGC 4826. We interpret the line ratios through ra- 
diative transfer modeling and map the physical conditions of the 
molecular gas. The main results of this study can be summarized 
as follows: 

- In both galaxies the three lines are similarly distributed over 
the data cubes, i.e. both spatially and in velocity. 

- In both galaxies the line ratios seem to be lower than in other 
galaxy centers. 

- In both galaxies the gas with greater than 0.2 is confined 
to Paa emitting regions detected by NICMOS. R21 instead 



reaches high values away from the Paa emitting regions de- 
tected by NICMOS. 

- According to a pseudo-LTE model, the molecular gas in 
NGC 4569 is cold, mainly optically thick in the CO(1-0) and 
(2-1) lines and with a fraction less than 50% optically thin 
in the CO(3-2) line. LVG modeling suggests the presence 
of an elongated ring of cold and dense gas coinciding with 
the ILR. There are two denser and colder clumps along this 
ring and one denser and warmer clump close (~160pc) to 
the nucleus. 

- More excited gas is resolved in the circumnuclear disk 
(CND) of NGC 4826. According to the pseudo-LTE model 
this corresponds to warmer gas with a ~50% of the CO(3-2) 
emission being optically thin. The LVG modeling indicates 
the presence of a semicircular arc of dense and cold gas cen- 
tered on the dynamical center and ~70pc in radius. The gas 
temperature increases and the density decreases toward the 
center. 

- NGC 4826 shows a near side/far side asymmetry in the CO, 
Paa and R32 distributions. This asymmetry suggests that the 
gas opacity could play a role in the CO morphology. It can 
be reproduced if the cold dense gas is distributed in a ring of 
70 pc radius embedded in a flared toroidal structure and can 
partly absorb the CO emission, or if the clumps forming the 
toroidal structure are inhomogeneously heated by the nuclear 
activity and are partly optically thick. 

- Several published CO maps of nearby active galaxies seem 
to show a slight CO emission excess on the far side of the 
nuclear region similar to the one observed in NGC 4826 sug- 
gesting that this could be a common property of active galax- 
ies. If confirmed, this would favor the interpretation based on 
the model with inhomogeneous clumps heated by the nuclear 
activity. 

This study demonstrates the feasibility of mapping the molecular 
gas physical properties in the centers of nearby galaxies based on 
multitransition CO observations with current interferometers. It 
also shows how important these observations are for correctly in- 
terpreting the CO emission and thus understanding the processes 
at work in nuclear regions. While the results derived in this 
study are still mainly qualitative and limited by uncertainties and 
resolution, they suggest new areas of investigation. Analyzing 
multitransition interferometric CO observations for more nearby 
galaxies with various inclinations should provide a statistical ba- 
sis for investigating the unified AGN paradigm. More impor- 
tantly, with ALMA, such observations will be greatly improved 
in terms of sensitivity (about one order of magnitude), resolution 
(one order of magnitude) and imaging quality (more complete uv 
coverage). Furthermore, it will be possible to observe more CO 
transitions with a single array and thus, to better constrain the 
models and to derive more quantitative results on the ISM prop- 
erties close to galaxy centers. For example it will be possible to 
keep the CO abundance as a free parameter and thus resolve the 
possible AGN influence on this parameter. It will also be possi- 
ble to resolve the clumpy structure of the CNDs and investigate 
in more detail the role of pressure in these structures. 
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Fig. A.l. Radial profiles of the molecular gas density (blue) and 
temperature (red) used in the models. 



Appendix A: Large Torus Models 

This appendix describes in more detail the two models discussed 
in Section l6*T2l We emphasize that these models are mainly geo- 
metrical toy models and they do not rely on any realistic simula- 
tions of the ISM. Their physical parameters were not fitted to the 
observations. These models are only meant to test whether large 
scale or clump scale opacity effects in a generic toroidal struc- 
ture could be consistent with the near side/far side asymmetry 
observed in NGC 4826. 

A.1. Model with large scale opacity 

In this model the clumps are assumed to be homegeneous and 
the only parameter related to the dumpiness of the medium is 
the volume filling factor. The sizes and shapes of the individ- 
ual clumps are not taken into account. The density and tempera- 
ture radial profiles used are represented in Fig. IA.ll The density 
increases exponentially with radius from « m [ n = 10 25 crrT 3 at 
lOpc to « max = 10 4 ' 9 crrr 3 at 70 pc and it decreases exponen- 
tially with height above the galactic plane. Inside 70 pc the tem- 
perature is related to the density by T oc n" with T m - dx — 50 K 
at lOpc and r m ; n = 10K at 70pc, i.e. a = -0.3. Outside 
70 pc the density decreases exponentially with radius down to 
"out = 10 2 cnT 3 and the temperature is kept constant at 12 K. 
Overall the density and temperature profiles define a cold dense 
ring at 70 pc. The gas volume filling factor is proportional to the 
gas density (where the gas is denser, it fills more of a given vol- 
ume). We assume that gas with a density higher than 0.6 x « max 
has a volume filling factor of 20% and can absorb the CO line 
emission from gas further away on the line of sight having the 
same line of sight velocity. We assume circular rotation and a 
gaussian line profile with a velocity dispersion of 15 km s out- 
side 70 pc and 30km s inside. The CO line intensities are com- 
puted with the LVG model. The resulting integrated CO maps 
and line ratio maps are presented in Fig. [23] 

A.2. Model with clump scale opacity 

In this model the individual clumps are assumed to be spheri- 
cal and inhomogeneous and they cannot obscure each other. In 
each clump the hemisphere exposed to the nucleus is assumed 
to be warmer and less dense than the other hemisphere. In ad- 
dition, the warm hemisphere is assumed to be optically thin in 
the CO(3-2) line, while the other hemisphere is assumed to be 
optically thick. CO(1-0) and (2-1) lines are assumed to be opti- 
cally thick everywhere. The density and temperature profiles of 
the warm hemisphere of the clumps as a function of the distance 



to the nucleus and the height above the plane are the same as in 
the previous model. The temperature profile of the cold hemi- 
sphere follows the same exponential profile with the same lower 
value but the maximum value is fixed to 30 K instead of 50 K. 
Thus, the difference in temperature between the warm and cold 
hemispheres increases with the temperature and reaches 20 K at 
the maximum. The density of the gas in the clumps is related to 
the temperature as in the previous model and the number density 
of clumps is assumed proportional to this gas density. The line 
intensities are computed assuming LTE. The resulting CO maps 
as well as the ratio maps are presented in Fig. [24] 



